Role of final state interaction and of three-body force 
on the longitudinal response function of 4 He 
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We present an ab-initio calculation of the longitudinal electron scattering response function off 
4 He with two- and three-nucleon forces and compare to experimental data. The full four-body 
continuum dynamics is considered via the Lorentz integral transform method. The importance of 
the final state interaction is shown at various energies and momentum transfers q. The three-nucleon 
force reduces the quasi-elastic peak by 10% for q between 300 and 500 MeV/c. Its effect increases 
significantly at lower q, up to about 40% at g=100 MeV/c. At very low q, however, data are missing. 
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Inelastic electron scattering off nuclei provides impor- 
tant informations on nuclear dynamics. Varying the mo- 
mentum q, transferred by the electron to the nucleus, one 
can focus on different dynamical regimes. At lower q the 
collective behavior of nucleons is studied. As q increases 
one probes properties of the single nucleon in the nuclear 
medium and its correlations to other nucleons from long- 
to short-range. Thus the inclusive longitudinal Rl and 
transverse Rt response functions are of particular impor- 
tance. Different from Rt, in a non-rclativistic framework 
Rl does not require the knowledge of implicit degrees of 
freedom (exchange currents), providing a clean leptonic 
probe of the nuclear Hamiltonian. In addition, the the- 
oretical study of inclusive processes is important to help 
planning further investigations, for selected kinematics, 
via exclusive scattering experiments. 

In the '80 and '90's an intense experimental activity 
has been devoted to inclusive electron scattering, (e,e'), 
in the so called quasi-elastic (q.e.) regime, corresponding 
to g-values of several hundred MeV/c and energy trans- 
fers uj around the q.e. peak (uj ~ q 2 /2m). Here one 
can envisage that the electron has scattered elastically 
with a single nucleon of mass m. Various nuclear targets 
have been considered, from very light to heavy ones 
At these q one enters a very challenging regime, where 
nuclear and subnuclear degrees of freedom inter wine. A 
very alive debate has taken place about the interpreta- 
tion of those data. The two most discussed topics have 
been: (i) short-range correlations, i.e. the dynamical 
properties of nucleons at short distances; (u) in medium 
modifications of the nucleon form factor. To date the 
debate is still open. More experiments are planned at 
Jefferson Laboratory (E05.110 at Hall A) which will con- 
tribute to those issues and a theoretical effort is needed 
to help interpreting old and new experimental results. 

The reason for concentrating on the q.e. regime has 
been the conviction that for such a kinematics the plane 
wave impulse approximation (PWIA) might be a reliable 
framework to describe the reaction. The neglect of the 
final state interaction (FSI) has the advantage to allow 



a simple interpretation of the cross section in terms of 
the dynamical properties of the nucleons in the ground 
state. Thus it is important to clarify the reliability of 
the PWIA (as well as of further refinements). The Eu- 
clidean approach has already shown that the PWIA is 
rather poor, however, this method does not easily allow 
to obtain the w-dependence of the FSI effects. 

The aim of this letter is twofold. On the one hand we 
study the role of FSI on R L of 4 He at 300 MeV/c < q < 
500 MeV/c, where by now only calculations with central 
two- nucleon forces exist [J H|. Here we use a realistic 
two-body potential augmented by a three-nucleon force 
(3NF) and compare the PWIA to results obtained via the 
Lorentz integral transform (LIT) method 0, Q . The LIT 
method is an ab-initio approach, which allows the full 
treatment of the four-body problem. It has already been 
applied to various realistic calculations of elcctroweak re- 
actions in three- 

@, 0, S 0, El and four-body systems 
f]~3l| . Different from the Euclidean approach, the LIT 
method allows a comparison with the PWIA regarding 
the w-dependence of Rl- Our second focus lies on the 
study of the role of 3NFs. We contribute to this much 
debated issue investigating 3NF effects on initial and fi- 
nal states by studying Rt, in various kincmatical regions. 

The choice of 4 Hc as a target is of particular interest. 
In fact 4 He has quite a large average density. Moreover 
its binding energy per particle is similar to that of heavier 
systems. Therefore 4 He results can serve better as guide- 
lines for investigating heavier nuclei than results for two- 
and three-body systems. Various inclusive 4 He (e. e') ex- 
periments have been performed in the past (see [14| for 
a summary of the world data), and a comparison theory- 
experiment is possible without the ambiguities, created 
by the Coulomb distortions, which affect heavier systems. 

The longitudinal response function is given by 

R L (u,q) = £ I <*/|p(<7)l*o> \ 2 S (ej + £j - E Q - J^, 
where M is the target mass, I'J'o//) an d ^0// denote ini- 
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tial and final state wave functions and energies, respec- 
tively. The charge density operator p is defined as 



r?) exp [iq • r; 



(1) 



where e is the proton charge and rf the isospin third 
component of nuclcon i. The <5-function ensures energy 
conservation. Rl contains a sum over all possible final 
states, which are excited by the electromagnetic probe, 
including also continuum states. Thus, in a straightfor- 
ward evaluation one would need to calculate both bound 
and continuum states. The latter constitute the major 
obstacle for many-body systems if one wants to treat the 
nuclear interaction rigorously. In the LIT method [H, Q 
this difficulty is circumvented by considering instead of 
R,L,(ui,q) an integral transform Cl (<j, q) with a Lorentzian 
kernel defined for a complex parameter a = /jr by 



£lO, q) 
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The parameter 07 determines the resolution of Cl and is 
kept at a constant finite value (07 ^ 0). The basic idea 
of considering Cl lies in the fact that it can be evaluated 
from the norm of a function which is the unique 

solution of the inhomogeneous equation 



(H-E -a)\^ q )=p(q)\* ) 



(3) 



Here H denotes the nuclear Hamiltonian. Due to the 
presence of the imaginary part 07 in (|3|) and the fact 
that its right-hand side is localized, one has a bound- 
state like asymptotic boundary condition. Thus, one can 
apply bound-state techniques for its solution. Finally, 
RL{oJ,q = const) is obtained by inverting the LIT (|2j). 
Subsequently the isoscalar and isovector parts of Rl are 
multiplied by the proper nucleon form factors. For the 
LIT inversion various methods have been devised (l5l.[l6j. 

The PWIA result is obtained under the hypothesis of 
one outgoing free proton with mass m and a spectator 
(A-l)-system with mass M s : 



dpn(p) S I u> 



(p + q) 
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Here n(p) represents the proton momentum distribution 
and e the proton separation energy. In the following we 
present results obtained with the Argonne V18 (AV18) 
[17t and the Urbana IX (UIX) [U] two- and three-body 
forces. As nuclcon form factor we use the proton dipole 
fit and the neutron electric form factor from [l9[ . The so- 
lution of ((3|), as well as the ground state \^o), is expanded 
in hypcrsphcrical harmonics (HH). The HH expansion is 
truncated beyond a maximum value K max of the HH 
grand-angular momentum quantum number. The con- 
vergence of the HH expansion is improved by intro- 
ducing a if maa; -dependent effective interaction (EIHH- 
method) [10, Hl[. In order to evaluate Cl we have cal- 
culated the norm ) directly, using the Lanczos 
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FIG. 1: Isovector (a) and isoscalar (b) parts of R.L(u>,q) at 
q=500 MeV/c. Single multipole contributions with K^ ax 
(solid) and with K^ x — 2 (dashed) first inverted and then 
summed up; single multipole contributions with K^ x first 
summed up and then inverted (dotted). 



algorithm [22j. The operator p is expanded in Coulomb 
multipolcs of order J. The LIT is calculated for each 
isoscalar (T=0) and isovector (T—l) multipole separately 
up to a maximal value of J ma x where convergence of the 
expansion is reached. The values of J m ax vary from 2 to 
7 for q ranging from 50 to 500 MeV/c. 

The accuracy of the results is determined mainly by 
the convergence of the HH expansion and the stability 
of the inversion. In the calculations we used a ground 
state hyperspherical momentum value K^ ax =16 (14) 
for the AV18+UIX (AV18) case, leading to a binding en- 
ergy of 28.4 (24.3) MeV. Since a multipole dependent 
convergence pattern has been encountered and each mul- 
tipole contributes differently to the total strength, the 
Kmax used for the LIT evaluation vary according to the 
value of J, namely K^ ax = 12 — 16 for even J and 
^max = 13 — 17 for the odd J have been considered. 



JT 
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Our LIT results converge at a percentage level. In Fig. [1] 
the accuracy of the results for Rl regarding both the 
HH expansion and the inversion stability aspects is il- 
lustrated exemplary for the isoscalar and isovector parts 
at g=500 MeV/c. The figures contain three curves: the 
full line is obtained when the single multipole contribu- 
tions C J L T ', calculated up to K^ ax , are first inverted and 
then summed up. The dashed line represents the results 
where the various multipole contributions C J L T are cal- 
culated only up to K^ ax — 2. The comparison between 
these two results illustrates the quality of the HH con- 
vergence. The dotted line reflects the inversion of the 
total Cl(&, q), where the various multipole contributions 
C'[ T , calculated up to Kf r ^ ix , are first summed up and 
then inverted. The comparison between the dotted and 
full lines shows the accuracy of the inversion. In Fig. 1 
one finds very satisfying results for both isospin channels 
for the HH convergence and the accuracy of the inver- 
sion as well. We should mention that we do not show the 
low-energy isoscalar response, where a narrow + reso- 
nance with a width of a few hundred keV is present at 
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FIG. 2: Rl(oj, q) at various q: PWIA using n(p) of 
AV18+UIX [2| (dotted); full calculation with AV18 (dashed) 
and AV18+UIX (solid). Data from Bates [25J (squares), 
Saclay [2(| (circles) and world-data set from [14| (triangles). 
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3: Rl(u,q) at various q with the AV18 (dashed), 



AV18+UIX (solid) and MTI-III (dash-dotted) potential. 
Data in (a) from [291 ]. 



TABLE I: Rl peak position ui p and Rl peak height without 
3NF (AV18), with 3NF (AV18+UIX), and relative 3NF effect 
AR = 100 x (i? L (AV18) - ,Rl(AV18 + UIX))/i? L (AV18). 
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E r very close to threshold [23|. To get accurate results 
for such a resonance a convergent LIT calculation with a 
07 much smaller than the presently used values (smallest 
value a i = 5 MeV) should be carried out, which then 
leads to a very slow asymptotically fall off of the solution 
I^Sq) (see [13]). Such a calculation requires a consider- 
able additional computational effort and thus the thresh- 
old region is excluded from our present work. Allowing a 
narrow resonance in the inversion [2?| . we have checked 
that our results are stable for energies above E r + 2ai. 

In Fig.[2]the results of Rl(w, q) at various q are shown 
and compared to data. In all cases one finds that the FSI 



effects are very large and essential for reaching agreement 
with experiment. The PWIA fails particularly in the q.e. 
peak and at low us. With growing q FSI effects decrease in 
the peak region, but not at low to. One may also consider 
a more refined PWIA, where a spectral function is used 
instead of a momentum distribution (see e.g. [28j|). In [28j 
it was shown that such an improved PWIA modifies the 
simple PWIA result by only 10-20 %. 

In Fig. [5] one also sees the 3NF effects on the full cal- 
culation. For <j=300 MeV/c one notes a good agreement 
of the data with the AV18+UIX result. This is true for 
g=400 MeV/c as well, if one does not consider the data 
of [25|, which exhibit larger error bars. At <7=500 MeV/c 
some discrepancies between theory and experiment are 
present in the low- and high-energy range, while there 
is a fairly good agreement in the peak region. However, 
investigations on the three-body systems [Kj have shown 
that a consideration of relativistic effects becomes impor- 
tant at such a momentum transfer. 

Table Q] illustrates the 3NF effect on peak position and 
peak height also for lower q. One notes that there is no 
unique 3NF effect on the position, while one has a reduc- 
tion of the height due to the 3NF at any q. The size of the 
reduction amounts to 10% for the higher q, whereas be- 
low g=300 MeV/c the reduction grows with decreasing q, 
reaching almost 30% at q <100 MeV/c. In Fig. |3] the re- 
sults at lower q are shown. The important role of the 3NF 
is evident in the whole peak region, leading to a strong 
decrease of Rl of up to 40% for some oj values. Recently 
also some new data at q ~ 200 MeV/c have been pub- 
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lished [29[ (see Fig. [3^,). While one finds a satisfactory 
agreement between the AV18+UIX result and data be- 
yond the peak, one observes a non negligible discrepancy 
in the peak itself. In Figs. [3^, [3]d we also illustrate Rl 
for a calculation |4j with a central two-nucleon potential 
(MTI/III model [301). Results are more similar to the 
AV18 than to the AV18+UIX curves, showing that the 
3NF effect is not simply explained by the binding energy 
difference ( 4 He binding energy with AV18, AV18+UIX, 
and MTI/III is 24.3, 28.4 and 30.6 MeV, respectively). 

We summarize our results as follows. We have carried 
out an ab-initio calculation of the longitudinal (e, e') re- 
sponse function RL(uj,q) of 4 He for various kinematics 
up to 5=500 MeV/c. The full dynamics of the four-body 
system has been taken into account for the 4 He ground 
state and the four-body continuum states as well. The 
rigorous inclusion of FSI has been achieved by use of the 
LIT method. Our work has been mainly focused on two 



points, namely the study of the importance of FSI and of 
3NF. We have shown that both ingredients play an im- 
portant role and need to be considered in a calculation 
of Rl- A particularly important finding are the very 
large 3NF effects of up to 40% in the Rl peak region 
at q < 200 MeV/c. Thus it is becoming apparent that 
there exists an electromagnetic observable, complemen- 
tary to the purely hadronic ones, where one can learn 
more about the not yet well established 3NF. In view of 
our findings we hope for a revival of the experimental 
interest in electron scattering, especially on light nuclei 
and at lower energies and momenta. 
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